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ABSTRACT 

It is now reasonably clear that disc fragmentation can only operate in the outer parts of proto- 
stellar discs (r > 50 au). It is also expected that any object that forms via disc fragmentation 
will have an initial mass greater than that of Jupiter. However, whether or not such a process 
actually operates, or can play a significant role in the formation of planetary-mass objects, is 
still unclear. We do have a few examples of directly imaged objects that may have formed in 
this way, but we have yet to constrain how often disc fragmentation may actually form such 
objects. What we want to consider here is whether or not we can constrain the likely popula¬ 
tion of planetary-mass objects formed via disc fragmentation by considering how a population 
of objects at large radii (a > 50) au - if they do exist - would evolve under perturbations from 
more distant stellar companions. We find that there is a specific region of parameter space to 
which such objects would be scattered and show that the known exoplanets in that region have 
properties more consistent with that of the bulk exoplanet population, than with having been 
formed via disc fragmentation at large radii. Along with the scarcity of directly-imaged ob¬ 
jects at large radii, our results provide a similar, but independent, constraint on the frequency 
of objects formed via disc fragmentation. 
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1 INTRODUCTION 


The most widely accepted mecha nism for the formation of plan¬ 
ets is the core accretion model dPollack et all 1996b . In this model, 
dust grains grow rapidly to form kilometre-sized planetesimals that 
then coagulate to form a rocky core dSafronovll 1 9721) which, if suf¬ 
ficiently massive, may then accrete a gaseous envelo pe to from a 
gas giant planet (Lissaue rll 19931: [ P ollack et alill996i ). An alterna¬ 
tive suggestion dKuiperll 195 ll : iBossll 1998h is that gas giant planets 
may form via direct gravitational collapse in discs that are suffi- 
ciently massive so as to sustain a gravitational instability dToomrd 

11964 . 


The advantage of the latter mechanism is that it en¬ 
sures that gas giant planets can form prior to the dispersal 
of the gas disc, thought to typically occur within ~ 5 Myr 
dHaisch, Lada & Ladal 1200111 . Most of the evidence, however, 
favours the standard core accretion mechanism. For example, 
gas giant planets are preferentiall y found around me t al-rich stars 
dSantos, Israelian & Mayor! |2004 iFischer & Valentil 1 2004 . and 
there is an in dication of a signature of the snowline in the exoplanet 
distri bution dSchlaufman. Lin & Idal 120091 : iRice. Penny & Home! 
2013). Both of these suggest that the amount of solid material in 
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the disc, and the distribution of this material, influences planet for¬ 
mation, which would not be expected if disc fragmentation were a 
dominant formation mechanism. 

Additionally, it is now fairly clear that disc fragmentation is 
physically implausible in the inner regions of protostellar discs 
dRafikovil2005it . Fragmentation requires that the disc be both g rav¬ 
itationally unsta ble, and that it be able to cool rapidly dGammiel 
I200ll : IRice et alj|2003l ). The inner regions of protostellar discs are 
likely too optically t hick to cool s ufficiently fast for fragm enta- 
tion to be possible dClar 5 120091 : IRice & Armitagel l2009h . The 
outer regions (beyond ~ 50 au), however, may well have condi¬ 
tions suitable for fragmentation and it has, consequently, been sug- 
ested that this may explain so me of the directly imaged exoplanets 
Kratter & Murray-Clay 20ld ). such as those in the HR8799 sys¬ 
tem 1 Marois et al.ll2008l). 

It has also been suggested dNavakshinl20 1 Pali that planets may 
form beyond 50 au via disc fragmentation, and then spiral inwards, 
resulting in them orbiting much closer to their parent stars than 
where they formed. Also, if both the grain sedimentation timescale 
and migration rate are fast, the planet may produce a core and lose 
some of its o uter envelope through tidal interactions with the parent 
star dN avakshhill20 1 0b| ). As such, this process could turn what was 
originally a massive gaseous planet, into a lower mass Neptune-like 
planet, or even a terrestrial planet. 

Recent population synthesis calculations dForgan & Rice! 
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l2013ah . however, suggest that this is very unlikely, and that it is 
much more likely that planets forming at large radii remain mas¬ 
sive and remain with much la rger radii than is typical for the 
known exoplanet population dMarcv et al.l2008l) . These initial pop¬ 
ulation synthesis models also ignored subsequent mass accretion 
onto the planets and, hence, the masses are lower limits; we’d ex¬ 
pect planets formed via disc fragmentation to have masses higher 
than these models suggest. An initial study to also estimate how 
such a population would evolve through dynamical interactions 
suggests that some (~ 25%) woul d be ejected, with others scat - 
tered onto high-eccentricity orbits (Forgan, Parker & Ricej[2015|). 
However, this work did not consider the subsequent evolution of 
these high-eccentricity planets through tidal interactions with their 

parent stars._ 

In this work we expand on iForgan. Parker & Rice! d2015h by 
considering how a population of planetary-mass bodies forming 
beyond 50 au would evolve through dynamical interactions with 
an outer popu latio n of stel lar com panions, that drive Kozai-Lidov 
oscillations dKozaill963 : lLidovl 1962b . and then through tidal inter¬ 
actions with the host star. What we’re aiming to do is to establish 
the orbital properties of planetary-mass objects that originated at 
large radii, and to determine if a population of such objects exists 
within the known exoplanet population. In Section 2 we describe 
the models that we use. In Section 3 we discuss our results, and we 
draw our conclusions in Section 4. 


2 BASIC MODEL 


The goal is to study the evolution of a planetary system that is being 
perturbed by a third, stellar-mass body on an outer orbit. Specifi¬ 
cally, we aim to establish if we can identify a population, within the 
known exoplanet population, that could have originated with large 
initial semi-major axes (a > 50 au) and that has then been scattered 
onto orbits with semi-major axes inside ~ 5 au. 

To do this, we use the equations that describe the sec- 
ular evolution of a star-p lanet-star system, first presented by 
lEggleton & Kiseleva! 12001 ) and which have been extensively used 
to mo del the evolution of exoplanets that are perturbed by a third 
body dWu & Murravl2003l ; lFabrvckv & Tremainel2007l ). The equa¬ 
tions include perturbing accelerations from a third body, general 
relativistic apsidal precession, perturbing accelerations from stel¬ 
lar and planetary distortions due to tides and rotation, and tidal in¬ 
teractions between the planet and its host star. We do, however, 
typically ignore the perturbing acceleration due to the planetary 
distortion as it only becomes significant when the planet is very 
close to its parent star, and re duces the timestep signifi cantl y. Here 
we use the form presented bv lBarker & Qgilvia d2009h and [Barker] 
d201 lh which are regular at e = 0, and we e xpand the contri¬ 
bution due to the third bo dy to octupole order dNaoz et alJl201ll : 
iNaoz. Farr & Rasiol 2012b . We also include a simple stellar wind 
model i Kawalei 1988b that allows the planetary host star to spin 


down as it e volve s on the Main-Sequence. Details of the model can 
be found in lRicel d2015b . 


2.1 Initial conditions 

To understand how a population of planetary-mass bodies that form 
in the outer parts of protostellar discs via disc fragmentation, are 
perturbed by an outer stellar companion, we run a Monte Carlo- 
type simulation in which we randomly choose the outer planet to 


have a semi-major axis between 50 and 80 au, and to have an ec¬ 
centricity that is initially small (drawn from the positive side of a 
Gaussian distribution with a half-width of e = 0.025). We don’t 
actually have a good sense of the initial properties for such a pop¬ 
ulation of planets. We could have selected the semi-major axis 
randomly in log a, but that the range is quite small (50 — 80 au) 
probably means that this wouldn’t make much difference, which is 
confirmed by a simple check of^our results. Also, the only simu¬ 
lations to consider this dStamatellos & Whitworth) [2009b suggests 
that randomly in a may be^a mor e reasona ble choice. Similarly, 
other simulations flHall. Forean & Ricell2015h suggest that the ec¬ 
centricity of those objects that survive at large radii is typically 
small. We also don't know the mass distribution for such a pop¬ 
ulation of planets, but since we’re interested in the possibility of 
planetary-mass companions forming at these radii, we assume the 
mass distribution is the same as th at of t he known population of 
exoplanets ( dN/dM p oc Mp 1 ' 15 ) dMarcv et al.ll2008h . The planet 
host star is assumed to have a mass of M* = 1 Mg. 

That 30 — 40% of Solar-like stars have stellar, or sub-stella r, 
companions dDuquennov & Mavoil 1 199 lb iRaghavan et alj 12010b . 
makes it likely that such companions will influence the evolution 
of the planetary system. We, therefore, assume that there is a com¬ 
panion with a mass of M 0 = 0.5 Mg and with a semi-major axis, 
chosen randomly in log a, b etween a 0 = 200 and a 0 = 20000 au 
dDuquennov & Mavon[l99lb . Both radial velocity and transit exo¬ 
planet searches tend to remove binary systems from their catalogs. 
However, these tend to be systems where the separation is smaller 
than we’ve assumed here dBrown et al]l2011a] ; Valenti & F ische r! 
1 2005b and so this should not bias our results. 

The companion’s eccentricity is chosen randomly to be be- 
twe en e 0 = 0 and e p = 1, although we impose stability crite- 


riaj Mardling & Aarsethll200il : lLithwick & Naoj|201 ll : Na oz et all 

120131) that ensure that the triple system is long-term stable and that 

the quadru pole and oc tupole terms in the secular equations dom¬ 
inate (see iRicel 120151 for details). We could choose to vary the 
mass of this outer companion, b ut since we are in the test-particle 
regime dLithwick & Naozll201 ll ). this should not influence the re¬ 
sults. We also fix the outer companion’s orbit to be in the xy plane 
and orient the inner orbi t so that the m ut ual inclina tion, i, is dis¬ 
tributed isotropically dWu. Murray & Ramshail 2007b . We also ran¬ 
domly orient the longitude of the planet’s ascending node, and the 
argument of periastron. 

Since we’re also interested in how the planet will evolve 
tidally with its host star (if scattered into an orbit that with 
a sufficiently small periastron) we assume that the star and 
planet have tidal quality factors of Q' = 5 x 10 6 and 

Q'„ = 5 x 10 5 iGoldreich & Soteill 19661; [Yoder & PealelfT98lL 


Jacfaqnj^jreenbem&jjunesl l2008t iBaraffe. Chabrier & Barmanl 
20ld ; Brown et alfeoilbl) . and tidal love numbers of k a = 0.028 


and k r = 0.51 dPetrovichll2015b respectively. W e also ass um e that 
any planet that reaches its Roche limit dFaber. Rasio & Willemsl 
120051) . given by. 


R„ 


M, 


0.462 \M, 


1/3 


( 1 ) 


is tidally destroyed. 


3 RESULTS 

As mentioned above, the goal is to establish if we can identify a 
population of exoplanets, from within the known population, that 
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Figure 1. The final semi-major axis distribution for a population of 5000 
planets, initially located between 50 and 80 au and then perturbed by an 
isotropically distributed population of stellar companions with semi-major 
axes between 200 and 20000 au. Most remain at large radii, but ~ 5% 
are perturbed into orbits with a < 1 au. A further ~ 17% are perturbed 
into orbits with periastra inside their Roche limit and are, hence, tidally 
destroyed. 


Figure 2. Figure showing the final eccentricity, plotted against semi-major 
axis, for those planets that are scattered onto orbits - and survive - with 
a < 5 au. The diamonds show our simulated planets, while the filled cir¬ 
cles are the known exoplanets with masses above 1 Jupiter mass, or radius 
above 1 Jupiter radius, and with a < 5 au. The planets scattered from large 
radii either become ‘hot’ Jupiters, with e ~ 0, or proto-hot Jupiters. This 
constrains the region of e — a space where such planets could be found. 


may have originated at large radii (a > 50 au) and that has then 
been scattered onto closer orbits via perturbations from an outer 
stellar companion that drives Kozai-Lidov oscillations. For our ini¬ 
tial simulations, we randomly select - as described above - the prop¬ 
erties of the inner and outer orbit and evolve the system, using a 
fourth-order Runge-Kutta integrator, for a randomly selected time 
of between 200 Myr and 10 Gyr. We then repeat this to produce a 
sample of 5000 systems. 

Figure [I] shows the final semi-major axis distribution and il¬ 
lustrates how the perturbations from the outer stellar companion 
sculpts the initial semi-major axis distribution. A large fraction - 
with a distribution that extends above the limits of the y-axis - re¬ 
main with semi-major axes between 50 and 80 au. Another pop¬ 
ulation undergoes Kozai-Lidov oscillations which produces highly 
eccentric orbits that lead to these planets tidally interacting with 
their parent stars, and ending up on orbits with semi-major axes 
a < 1 au. Another population (no t sho wn in FigureQ} reaches their 
Roche limit (iFaber, Rasio & Willems!IT005i ) and is assumed to be 
tidally destroyed. 

From an initial population of 5000 systems, our simula¬ 
tions suggest that 864 (~ 17%) are tidally destroyed, 237 (~ 
5%) survive inside 1 au, and the rest primarily remain out be¬ 
yond 50 au. Binary companions to Solar-type stars are quite 
com mon, with 30-40% having stellar, or sub-stella r, compan¬ 
ions touquennov & Mavotlll99l1 : lRaghavan et al.l2010h with semi¬ 
major axes between 20 and 20000 au. We’ve assumed that these 
companions have semi-major axes that are evenly distributed in log 
a and have only considered companions with a > 200 au. This 
would suggest that 20 —30 % of Solar-type stars might satisfy these 
condition. Given that these results suggest that ~ 5 % of such sys¬ 
tems could scatter an outer planet (a initially between 50 and 80 
au) into an orbit that gets tidally circularised - and survives - inside 
a = 1 au, would then suggest that 1 - 1.5 % of all Solar-type stars 
could have such close-in planets if all such stars form planetary- 
mass companions beyond 50 au via disc fragmentation. 

Admittedly, we’ve only presented results from a single set of 
parameters here. Flowever, we have run some tests with different 
Q' p and Q' s values and the results are broadly similar. For very large 


values of Q' p ( Q' p > 10 7 ) the fraction surviving inside 1 AU can 
drop below 4 %, but that would still suggest that ~ 1 % could have 
such close-in planets if companions that form at large radii via disc 
fragmentation are very common. 

Figure [2] shows the eccentricity, plotted against semi-major 
axis, for all those planets that have final semi-major axes inside 
a = 5 au (diamonds) and also shows (solid circles) all the known 
exoplanets with masses above M p = 1 Mj up or - for those with¬ 
out a mass estimate - with radii above R p = 1 Rj up . What it 
shows is that those planets scattered from beyond 50 au will al¬ 
ways end up either as ‘hot’ Jupiters wit h a < 0.1 au and with sma ll 
eccentricities, or as proto-hot Jupiters jDawson & Johnsonil2012) ). 
Proto-Hot Jupiters are those gas-giant planets that are found on the 
high-eccentricity boundary, beyond which planets would typically 
be tidally destroyed, and that will ultimately tidally evolve to be¬ 
come ‘hot’ Jupiters. This is quite useful in that it constrains the 
region of e — a space where we would find such planets, if they do 
indeed exist. 

Again, we’ve only presented results from a single set of pa¬ 
rameters. We’ve used a relatively low Q p value and we find that 
higher Q' p values tend to result in proto-hot Jupiters having larger 
eccentricities, for a given final semi-major axis. In fact, for suffi¬ 
ciently large Q p values ( Q p > 10 7 ) the proto-hot Jupiters end up 
in a region of e — a space where very few exoplanets are found. It 
does seem clear, therefore, that if planetary-mass bodies are scat¬ 
tered from very large initial radii, they would likely be found either 
as ‘hot’ Jupiters or as proto-hot Jupiters, near the boundary beyond 
which they would probably be tidally destroyed. 

Figure [2] introduces the first possible issue with this mech¬ 
anism playing a significant role in producing close-in exoplan¬ 
ets. As mentioned above, if planet formation via disc frag¬ 
mentation is common (i.e., most stars have gas giants orbit¬ 
ing beyond 50 au) then 1 — 1.5 % of stars could have close- 
in exoplanets that were scattered in from these large initial 
radii. However, Figure [2] shows that these would have to be 
either ‘hot’ Jupiters, or be proto-hot Jupiters undergoing high- 
eccentricity tidal migrati on. Currently, only about 1 % of Sun-like 
stars have such planets jFressin et al. 20131 : IPetigura et al.[|2013l : 
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IDawson, Murrav-Clav & Johnsonll2015t) and, so, if planet forma¬ 
tion via disc fragmentation is common, tins scattering scenario 
could explain almost all the known ‘hot’ and proto-hot Jupiters. 

However, one can also largely explain the ‘hot’ Jupiters and 
proto-hot Jupiters through a combination of disc migration and 
scattering (both from other planets and from stellar and sub-stellar 
companions) of planets that initially form at modest semi-major 
axes (inside ~ 10 au) jNagasawa, Ida & Besshq 120081) . In fact, 
we reran our simulations with the only change being that the 
planets were randomly distributed between 5 and 10 au, rather 
than between 50 and 80 au. In this case, just over 5 % were 
scattered onto close-in orbits, either ‘hot’ Jupiters or proto-hot 
Jupiters. This is possibly a slight under-estimate as we’ve re¬ 
stricted the companions to be beyond 200 au. H owever , i f 20 % 
of Sun-like stars have gas giant companions jMarcv et al.ll2008h . 
and 20 — 30 % of such systems have stellar companions, then this 
would suggest that 0.2 - 0.3 % of such systems will end up as ‘hot’ 
Jupite rs or pr oto-hot Jup iters vi a this mechanism alone. Moreover, 
when iDawson, Murrav-Clav & Johnsonl J2015h analyzed the dis¬ 
tribution of transit durations, and therefore eccentricities, among 
'hot' Jupiters found by Kepler , they found that eccentric proto-hot 
Jupiters are extremely rare compared to ‘hot’ Jupiters, indicating 
that the majority of the observed ‘hot" Jupiters likely arrive by disc 
migration and not by high eccentricity scattering from the snow¬ 
line or beyond. Therefore, that closer stellar companions, other 
planetary companions, and disc migration could also play a role 
in forming such planets, suggests that such mechanisms are suffi¬ 
cient to explain the observed population, and that scattering planets 
from beyond 50 au is unlikely to play a significant role. 

3.1 The properties of close-in exoplanets 

One way to further investigate this is to consider the properties of 
the known exoplanets that lie within the region to which these outer 
planets could be scattered. Figure [2] suggests that planets scattered 
from beyond 50 au can only become either ‘hot’ Jupiters, or proto- 
hot Jupiters. Consequently, we assume that only those known exo¬ 
planets lying to the left of the dashed line in Figure [2] could have 
been scattered from initial radii beyond 50 au. As mentioned above, 
the exact location 'me — a space to which proto-hot Jupiters could 
be scattered does depend somewhat on the tidal dissipation param¬ 
eters. However, using a higher Q' p value would likely move the 
boundary to slightly higher e values, and so our choice of bound¬ 
ary (dashed line in Figure [3} should at least capture most objects 
that could have been scattered from a large initial semi-major axis. 

Figure[3]shows the mass distribution for the known exoplanets 
with masses above 1 Jupiter mass that exist to the left of the dashed- 
line in Figure [2] (solid line) and also shows the mass distribution 
of all known exoplanets with masses above 1 Jupiter mass. Both 
distributions are normalised with respect to their maximum value. 

The two distributions are very similar and is therefore consis¬ 
tent with most ‘hot’ Jupiters and proto-hot Jupiters coming from the 
main exoplanet population. Additionally, most of the planets have 
masses below 5 Jupiter masses, and more than 65 % have masses 
b elow 3 Jupiter masses. It has been suggested dForean & Rice! 
1201 ih that the initial mass of planets that form via disc fragmenta¬ 
tion will typically exceed 3 Jupiter masses. Even after tidal down¬ 
sizing, th e vast majorit y of objects still have masses above 5 Jupiter 
masses dForean & Ricell2013al) and, i f there is any external irradia¬ 
tion, the initial mass increases dForean & Ricell2013bll . However, a 
large fraction of the observed ‘hot’ Jupiters and proto-hot Jupiters 
have masses below 3 Jupiter masses and, therefore, have masses 
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Figure 3. The solid-line shows the mass distribution of all known exoplan¬ 
ets with masses above 1 Jupiter that exist in the region to which outer 
planets could be scattered (found to left of the dashed-line in Figure [2}- 
The dashed-line shows the mass distribution of all known exoplanets with 
masses above 1 Jupiter mass. It’s clear that the two distributions are very 
similar and is consistent with ‘hot’ Jupiters and proto-hot Jupiters coming 
primarily from the bulk exoplanet population. 



below that expected for planets that formed via gravitational insta¬ 
bility. 

As discussed earlier, our simulations assumed that the objects 
forming at large initial radii had the same mass distribution as the 
bulk exoplanet population. Since it seems likely that objects form¬ 
ing via disc fragmentation would tend to have initial masses of 
a few Jupiter masses, this may not be an appropriate mass distri¬ 
bution. However, even if we repeat our simulations with the mass 
fixed at 15 Mj up , the fraction that survive inside 1 au drops from 
~ 5 % to ~ 2.5 %, and the region of e — a space to which they 
are scattered remains unchanged. Therefore, even if our mass dis¬ 
tribution is not quite representative of what would be expected for 
objects forming at large radii via disc fragmentation, this wouldn’t 
change our results significantly. 

It is also generally recognised that the host stars of gas gi¬ 
ants ten d to be metal-rich compared to a typical sample of simi - 
lar stars l iSantos, Israelian & Mavoj2004l ; lFischer & Valent J2005I) . 
This is thought to indicate that gas giant planets are more likely to 
form in discs that are enhanced in solids and i s regarded as c onsis- 
tent with the standard core accretion scenario (1 Pol lack et al.ll 19961) . 
Disc fragmentation would not seem to require this enhancement 
and, if anyt hing, it may even be more effectiv e in systems that are 
metal poor dMeru & Bate 2010; Clarke 200^) since such discs can 
cool more efficiently than discs that are enhanced in solids. It has 
been suggested, however, that there may be a similar metallicity 
dependence for planets that form via gravitational instability, un- 
dergo pebble accr etion, and then migrate rapidly to the inner disc 
dNavakshinl 120151) . This, however, does not necessarily apply to 
those that still have large semi-major axes once the disc has dis¬ 
sipated, and it is this population that is relevant to what we’re pre¬ 
senting here. 

Figure [4] shows the metallicity distribution for those known 
exoplanets in the region to which outer planets can be scattered 
(solid line) compared to all known exoplanets (dashed line). In both 
cases, we restrict the mass to be above M = lMj up or - for those 
without a mass estimate - the radius to be above R = 1 Rj up . It’s 
clear that the distributions are very similar and that - like the bulk 
exoplanet population - ‘hot’ Jupiters and proto-hot Jupiters tend 
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Figure 4. The solid-line shows the metallicity distribution for all known ex¬ 
oplanets - with masses above 1 Jupiter mass, or radii above 1 Jupiter radius 
- that exist in the region to which outer planets could be scattered (to the left 
of the dashed-line in Figure[2] The dashed-line is the metallicity distribution 
for all known exoplanets with masses, or radii, greater than that of Jupiter. 
That both populations are preferentially found around metal-rich stars is 
consistent with 'hot’ Jupiters and proto-hot Jupiters coming primarily from 
the bulk exoplanet population. 


to be found around metal-rich stars. Again, this is consistent with 
such planets being scattered, or migrated, from the bulk exoplanet 
population, rather than from radii beyond 50 au. 

To further quantify if the ‘hot’ and proto-hot Jupiters could 
have originated from the bulk exoplanet population, we carried out 
Kolmogorov-Smirnoff tests on the metallicity and mass distribu¬ 
tions. A Kolmogorov-Smirnoff test of the metallicity distribution 
returns a probability of Pks = 0.78, indicating that we can’t rule 
out that ‘hot’ and proto-hot Jupiters are drawn from the same distri¬ 
bution as the bulk exoplanet population (dashed-line in Figure [4). 
Similarly, the KS test of the mass distributions (Figure[3j returns a 
probability of Pks = 0.22, again indicating that we can’t rule out 
that the two populations are drawn from the same distribution. 

Figure[3]therefore suggests that most ‘hot’ Jupiters and proto- 
hot Jupiters have lower masses than would be expected for plan¬ 
ets forming via disc fragmentation, and the KS test indicates that 
we can’t rule out that they’re drawn from the same distribution 
as the bulk exoplanet population. This would appear inconsistent 
with a significant fraction originating from an outer population that 
formed via disc fragmentation. That we also can’t rule out that the 
metallicity distribution is drawn from the same distribution as the 
bulk exoplanet poopulation would also suggest that ‘hot’ and proto- 
hot Jupiters are more likely to have come from the main exoplanet 
population, than from a population of planetary-mass bodies that 
formed via disc fragmentation at large orbital radii. 


3.2 KELT1-B - The closest brown dwarf companion 

Possibly one of the more interesting obje cts in the c ontext of 
what we’re considering here is KELT-IB dSiverd et all 1201 A . It 
is a brown dwarf, rather than a planetary-mass, companion to 
a 1.335 Solar-mass star on an orbit with a semi-major axis of 
a — 0.024 au, and an eccentricity of e = 0.001. It has a mass 
of 27.4 Jupiter masses and so is more consistent with what might 
be e xpected for objects that form via disc fragmentation a t large 
radii dStamatellos & Whitworthl200^ : lForgan & Ricel2013al lbh. We 
didn’t directly address this system, as it is beyond what we can con¬ 


figure 5. Figure showing the final eccentricity, plotted against semi-major 
axis, for objects with a mass of M = 30 Jupiter masses, that started with 
a semi-major axis between a = 50 and a = 80 au, and that were scat¬ 
tered by an even more distant stellar companion. As with the simulations 
of planetary-mass objects, the only objects that survive inside a = 10 au 
are those that end up inside a = 0.1 au and have a small eccentricity, or 
those that are on the high-eccentricity boundary, beyond which they would 
be tidally destroyed. The location of KELT-1B is indicated at a = 0.024 
au and e = 0.001, showing that scattering from a large initial radius could 
potentially form such a system. 


sider here, but we did repeat our simulations with the only changes 
being that we simulated 2500 systems (rather than 5000) and as¬ 
sumed that the outer non-stellar object had a mass of 30 Jupiter 


As with our earlier simulations, a large fraction of the objects 
are scattered onto orbits that pass close to the parent star. In the 
simulations with an outer object with M = 30 Mj up , 558 (22 %) 
are tidally destroyed and 45 (1.8 %) survive on orbits with a < 1 
au. This means that a smaller percentage survive with small semi¬ 
major axis, than for the case where the outer object was assumed 
to be of planetary-mass, but a similar fraction are scattered. Again, 
as shown in Figure[5] they’re all either in the same location as ‘hot’ 
Jupiters or proto-hot Jupiters. beyond which they would be tidally 
destroyed. The position of KELT-IB is indicated at a = 0.024 au 
and e = 0.001. As is clear, it is possible for a system such as 
KELT-IB to form via scattering from a large initial radius. 

As already mentioned, abo ut 1 % of Sun-like stars have a ‘hot’ 
Jupiter o r a proto- hot Jupiter (F ressi n et al~l l20lH ; IPetigura et al.1 
l2013l : lDawson. Murrav-Clav & Johnsonl 2015t ). If we consider the 
known exoplanets with mass estimates, there are 127 with masses 
above 1 Jupiter mass and that are either ‘hot’ Jupiters or proto-hot 
Jupiters (i.e., lie to the left of the dashed line in Figure^. If we 
also consider those without mass estimates, but with radii above 
1 Jupiter radius, it increases to 282. KELT-IB is one of the few 
brown dwarf companions that have orbital properties similar to that 
o f a typic al ‘ hot’ J upiters. The only other known one is Corot-3b 
iTriaud et al.ll2009l) . a 21 Jupiter mass companion to a 1.36 Solar- 
mass star, with an orbital semi-major axis of a = 0.057 au. As 
mentioned above, a 30 Jupiter mass object, scattered from a large 
initial radius, is a few times less likely to survive inside a = 1 au, 
than a planetary-mass body. Consequently, if KELT-IB is indica¬ 
tive of a population of bodies at large orbital radius, a few of the 
known ‘hot’ and proto-hot Jupiters may have originated via scat¬ 
tering from such a population. In other words, maybe ~ 1% of the 
known ‘hot’ and proto-hot Jupiters originated from beyond ~ 50 
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au. However, as mentioned earlier, all the known ‘hot’ and proto- 
hot Jupiters could have such origins if just over half of Solar-like 
stars have a population of outer planetary companions. Therefore, if 
only ~ 1% of known ‘hot’ and proto-hot Jupiters have such an ori¬ 
gin, less than ~ 1% of Sun-like stars can have planetary, or brown 
dwarf, mass companions at large radii. 


3.3 Properties of wide exoplanets probed by direct imaging 


As we’ve shown above, planetary-mass bodies scattered from large 
initial radii (a > 50 au) can only be found as ‘hot’ Jupiters or proto- 
hot Jupiters. We’ve also shown that the properties of the known 
exoplanets in those regions are very similar to that of the bulk exo¬ 
planet population, suggesting that very few - if any - were scattered 
from beyond 50 au. A simplistic analysis based on the existence of 
a KELT-IB, a brown dwarf with orbital properties similar to that of 
'hot' Jupiters, does suggest, however, that maybe a small number 
of the known ‘hot’ and proto-hot Jupiters could have originated be¬ 
yond ~ 50 au. This would indicate that maybe ~ 1% of Solar-like 
stars have distant planetary, or brown dwarf, companions. 

Most planetary detection methods (e.g. radial velocity and 
transit, for instance) are, however, more sensitive to close-in 
planets (<5 AU) than to those that may be further out. In 
contrast, direct imaging is more sensitive to planets at wider 
separations, and of all planetary detection methods uniquely 
probes planets at large radii (>50 au). To date, ~20 wide ex¬ 
oplanet or very-low-mass objects (<25 Mj uv ) have been dis- 


-mass objects ( <2b M j uv ) have be e n dts- 
covered iMarois et al.l 120081 120101: lLaeranee et aD 12001 l20ld: 
iKuzuhara et alj 20131 among others). In the last decade, nu¬ 


merous 8m telescope + coronagraph surveys have searched for 


planetary-mass bodies around young stars (e.g. Lowrance et al. 

20051: Billet etal. 20071: 

Lafeniere et al. 20071; Chauvin et al. 

201CJ: Heinze et al. 201& 

Vigan et all 20121; Biller et al. 

2013|; 

Rameau et al. 20131: Wahhai et alj 20131; Chauvin et alj 

20 ii 

Brandt et al. 2014; Bowler et al. 2013), in order to detect new com- 


panions and also to determine the frequency of such wide compan¬ 
ions. Although individual surveys tend to only cover 80-100 stars, 
these surveys have - overall - now observed a sample of ~500 stars. 

They have unearthed few new planets and most have yielded 
null planet detections, leading to the unambiguous result that 
planetary-mass bodies are rare at large radii (>50 AU). For in¬ 
stance, Biller et al. (2013) find, from a Bayesian analysis of a 
null planet detection from 78 young moving group stars observed 
with the Gemini NICI planet-finder, that the frequency of 1-20 
Mj up companions at semi-major axes from 10-150 AU is <18% 
at a 95.4% confidence level using DUSTY models (Chabrier et al. 
2000) and is <6% at a 95.4% confidence level using COND mod¬ 
els (Baraffe et al. 2003). Brandt et al. 2014b combine imaging data 
from multiple surveys to build a 250 star survey with 5 detections 
of brown dwarf companions. Modelling this population with a sin¬ 
gle power-law distribution, they find that 1-3% of stars (68% confi¬ 
dence) host 5-70 Mj up companions between 10-100 AU. They ar¬ 
gue that this suggests that most wide companions formed through 
disc fragmentation rather than through core accretion. However, 
whichever formation mechanism is responsible for their formation, 
wide giant-exoplanet companions are quite rare. 

Further analysis, however, suggests that the few known wide 
companions are indeed more likely to have been formed by disc 
fragmentation, than via core accretion. Biller et al. 2013 marginal¬ 
ize over a wide range of potential planetary distributions and Brandt 
et al. 2014b explicitly fit power-laws unconnected to radial veloc¬ 
ity studies. However, most statistical treatments of directly imaged 


planetary companion distributions adopt fixed power-law distribu¬ 
tions in semi-major axis and planet mass drawn from radial ve¬ 
locity planet studies (e.g. Cumming et al. 2008) and extend these 
proscriptions to separations appropriate for directly imaged plan¬ 
ets (e.g. >10 AU). As the population of plan ets probe d by ra¬ 
dial v elocity is likely formed by core-accretion dFischer & Valentil 
120051) . these simple empirical power-laws have been interpreted 
(even when extended to wide separations) as the expected popu¬ 
lation as produced by core-accretion. Since the power-laws used 
in radial velocity studies extend to arbitrary semi-major axes, these 
studies add a new “cutoff’ parameter, truncating the power-law dis¬ 
tribution at a specific semi-major axis (assumed to be similar to the 
size of the primordial disc). These studies have found that the semi¬ 
major axis cutoff for core-accretion planets must be <75 AU (e.g. 
Nielsen and Close 2010) and likely less than 30-40 AU. 

Although this analysis suggests that planetary-mass bodies at 
large radii probably formed via disc fragmentation, rather than via 
core accretion, the direct imaging surveys are also consistent with 
the results of the simulations presented here: planetary-mass bodies 
at large radii are rare and, therefore, disc fragmentation rarely forms 
such objects. It is possible that a few percent of stars host such 
companions, but that is consistent with the suggestion in Section 
3.2, that at most a few of the known ‘hot’ Jupiters and proto-hot 
Jupiters could have been scattered from a large initial radius. 


4 CONCLUSIONS 

In this study we have aimed to investigate if we can identify a popu¬ 
lation, from within the known exoplanet population, that may have 
originated from large initial radii (a > 50 au). Given that a reason¬ 
able fraction (~ 30 — 40 %) of Solar-like stars have stellar compan¬ 
ions, a population of planetary-mass companions beyond 50 au - if 
it does exist - should be susceptible to scattering via perturbations 
from these even more distant stellar companions. The analysis here 
suggests that if any such planets are scattered to closer orbits they 
should be found as ‘hot’ Jupiters (a < 0.1 au and e ~ 0) or proto- 
hot Jupiters; gas-giant planets on the high-eccentricity boundary, 
beyond which they would typically be tidally destroyed. However, 
when we consider the known exoplanets that lie in that region of 
parameter space, they have properties (mass and metallicity distri¬ 
butions) that are consistent with that of the bulk exoplanet popula¬ 
tion. Similarly, the typical mass of the known exoplanets in these 
regions is less than that expected for planets that form via disc 
fragmentation at large radii. Additionally, direct imaging searches 
for planetary-mass bodies around Solar-like stars also suggest that 
such objects are rare. Therefore, we conclude that there is probably 
not a substantial population of planetary-mass bodies at large radii 
around Solar-like stars and, consequently, that disc fragmentation 
rarely forms such objects. 

Of course, there are a number of caveats to the above. We only 
presented results from a single set of parameters in our simulations, 
but most were chosen to be conservative. For example, the stellar 
wind parameter was chosen to spin the star down slightly more than 
may be reasonable. A weaker wind would remove less angular mo¬ 
mentum, slow the tidal evolution of the planet, and allow a larger 
fraction to survive. The tidal quality factors were also chosen to be 
on the low side of the range. Additional tests suggest that larger 
Q' p values wouldn’t change the fraction surviving inside 1 au sub¬ 
stantially and would tend to move the proto-hot Jupiters to even 
larger eccentricities than considered here. Therefore, if planetary- 
mass objects at large radii were common, they could explain a large 
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fraction of the known ‘hot’ Jupiters and proto-hot Jupiters. That 
such known exoplanets appear more consistent with an origin from 
within the main exoplanet population suggests that there is not a 
substantial population of planetary-mass bodies at large radii. 


As suggested above, a possible consequence of this analy¬ 
sis is that disc fragmentation rarely forms planetary-mass bod¬ 
ies. One possibility is that any object formed via this mecha¬ 
nism would quickly grow in mass to become a brown dwarf 
(M > 15 Mj up ) as sugg ested by disc fragmentation simulations 
Jstamatellos & Whitworth|[2009l) . An example might be KELT-1B, 
which has a mass of M = 27.4 Mj up and an orbital radius of 
a = 0.024 au. Our analysis suggests that it may indeed have been 
scattered from a large initial orbital radius. Additionally, KELT1-B 
is one of only a small number of known brown-dwarf mass objects 
that lie in the region of parameter space to which we’d expect outer 
objects to be scattered. Therefore, even if KELT-IB is an example 
of an object scattered from beyond 50 au, it would still be consis¬ 
tent with such objects being rare and, hence, that disc fragmentation 
rarely operates. 


Studies _ of _ s elf-gr a vitating _ disc _ ev olution 

dWhitworth & Stamatello j 120061 : iRice & Armitagel 120091 : Idarkel 
120091 : Forgan & Ricel 201 lh do, however, suggest that any disc 
with an initial outer radius that extends beyond ~ 50 au should be 
susceptible to fragmentation and, consequently, the formation of 
planetary or brown dwarf-mass objects. That disc fragmentation 
seems unlikely, may suggest that discs with such large initial 
radii are rare, or that external irradiation is ab le to stabilise such 
discs against fragmentation dRice et alj 1201 lh . There are indeed 
observations dMaurv et alJl2Qloh suggesting that discs around in 
very young Class 0 protostellar systems do tend to be compact, 
with outer radii inside 50 au. If so, then by the time such discs 
have expanded to larger radii, there may_be in suffi cien t mass for 
them to be susceptible to fragmentation dRice. Mayo & Armitagel 
l20ld) . There is, however, some evi dence for extend ed discs 
in very-young protostellar systems JTobineUtL 2012h. w hich 
may be gravitationally unstable dForgan & Ricel 12013d') . If 
such discs do exist, then the possibility that disc fragmentation 
rarely operates, may indicate that some thing, such as magnetic 
fields tCommercon. Henne belle & Henning 2011), is inhibiting 
fragmentation in such discs. 


Additionally, fragment destruction is a common out¬ 
come for all disc fragmentation studies, especially as self 
gravitating disc migration appears to be relatively rapid, 
and in duces tidal disruption bef ore the fragment is full y 
bound dB aruteau. Meru & Paardekooperl 1201 ll : IZhu et all 12012 : 
iForgan & Rice 2013ah . Dynamical evolution while still in the 
nascent stel lar cluster, could also remov e objects that survive at 
large radii dForgan. Parker & Ricel l20f^ ). This, however, would 
need to be very efficient if it were to reduce an initially large popu¬ 
lation of outer planetary-mass bodies, to one more consistent with 
the results presented here. It is possible, therefore, that the results 
here have implications for planet formation, suggesting that most 
known exoplanets formed via core accretion, and for star forma¬ 
tion itself, suggesting - for example - that discs around very young 
stars are typically compact, with initial outer radii inside the radius 
where disc fragmentation becomes viable, or that magnetic fields 
and external irradiation can stabilise radially extended discs against 
fragmentation. 
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